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O ! Abstract 

o ■ 

CN ' We study the exclusive decays of B meson into P-wave charmonium states Xcj(J 

0. 1) in the QCD factorization approach with light-cone distribution functions describing 
H-j ! the mesons in the processes. For B — ► XdK decay, we find that there are logarithmic 

£T) ' divergences arising from nonfactorizable spectator interactions even at twist- 2 order and 

CN . the decay rate is too small to accommodate the experimental data. For B — ► XcoK decay, 

\ we find that aside from the logarithmic divergences arising from spectator interactions 

' at leading-twist order, more importantly, the factorization will break down due to the 

infrared divergence arising from nonfactorizable vertex corrections, which is independent 
of the specific form of the light-cone distribution functions. Our results may indicate that 
QCD factorization in the present form may not be safely applied to 5-meson exclusive 
■ decays to charmonium states. 

o ; 
^ ■ 

Oh; PACS numbers: 13.25.Hw; 14.40.Gx 

& : 

Exclusive nonleptonic 5-meson decays provide a important opportunity to determine the 
^ ■ parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, to explore CP violation and to 
observe new physics effects. Recently, B physics has received extensive experimental attention 
such as from high-energy experiments at the Tevatron and at the e + e~ B factories. On the other 
hand, although the underlying weak decay of the b quark is simple, quantitative understanding 
of nonleptonic i?-meson decays is difficult due to the complicated strong-interaction effects. 

Beneke et al. have considered two-body nonleptonic 5-meson decays extensively including 
light-light as well as heavy-light final states within the QCD factorization approach^ El IS]. 
The general idea is that in the heavy quark limit ^> Aqcd ; the transition matrix elements 
of operators in the hadronic decay B —>■ MiM 2 with Mi being the recoiled meson and M 2 being 
the emitted meson can be calculated in an expansion in l/m^ and a s . The leading term in 
l/m& assumes a simple form[2j: 

(M 1 M 2 \O l \B)=F BAh (m 2 2 ) f duT I (u)<p M2 (u)+ ^dudv^^u^MOMv^iu), (1) 

Jo " Jo 

where M 2 is a light meson or a quarkonium and F BMl is the B — > Mi transition form factor; 
4>m is meson light-cone distribution amplitude and T 1 ' 11 are perturbatively calculable hard 
scattering kernels. If we neglect strong interaction corrections, eq.Q reproduces the result 
of naive factorization. However, hard gluon exchange between M 2 and BMi system implies a 
nontrivial convolution of hard scattering kernels T I,!I with the light-cone distribution amplitude 
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4>m 2 - This method works well for light-light final states^ El EJ El as well as heavy-light final 
states El- 

Exclusive 5-meson decays to charmonium are important since those decays e.g. B — > J/ipK 
are regarded as the golden channels for the study of CP violation in B decays. It is argued 
that because the size of the charmonium is small(~ l/a s m^) and its overlap with the (B,K) 
system may be negligible 0, the same QCD factorization method as for B — > tttt can be used 
for B — > J/ipK(K*) decay. Indeed, explicit calculations [7J E] show that the nonfactorizable 
vertex contribution is infrared safe and the spectator contribution is perturbatively calculable, 
where the light-cone distribution functions are used for B, K as well as J/ip mesons. This 
small size argument for the applicability of QCD factorization for charmonia is intuitive and 
interesting, but it needs verifying for charmonium states e.g. the P-wave Xcj states other than 
the J/i/j. In addition, recently BaBar and Belle collaborations have measured the exclusive 
decays of B — > Xc.jK(J = 0, 1)0 CD]- So, it is also interesting to compare the predictions based 
on the QCD factorization approach with the experimental data. In this letter, we report the 
problems of the QCD factorization approach encountered in these two decay channels. As in 
0E]> m the following we will use light-cone distribution functions to describe B, K, as well as 
charmonium mesons. 

We first consider B — > x&K decay. The effective Hamiltonian is written as[TT] 



p 10 

H eS = -| (v a V* (Cid + C 2 2 ) - V tb V* £ C& 

V 2 i=3 



where Gp is the Fermi constant, Cj are the Wilson coefficients and V qiq2 are the CKM matrix 
elements. Here the relevant operators Oi are given by 

Ox = {s a bp)v_ A ■ {cpc a ) v _ A , 2 = (s a b a ) v _ A ■ {cpC^v-A, 

03(5) = (s a b a )v-A ■ ^2(Qi3ll3)v-A(V+A), ^4(6) = (s a bp)v-A 4 ^2(Qfi<3a)v-A(V+A), (3) 

<? 9 

3 3 

^7(9) = -{s a b a ) V -A ■ X^^^V+AtV-A), 08(10) = ^(s<*Mv-A * e q(^l3 ( la)v+A(V-A), 
1 q 1 q 

where a, (3 are color indices and the sum over q runs over u,d,s,c and b. Here (qiq2)v±A = 
gi7 M (l ±7s)<?2. 

To calculate the decay amplitude, we introduce the Xci decay constant f Xel as[T2] 

(Xci|c(0)7 M 7 5 c(0)|0) = -if Xcl m Xcl e;, (4) 

where m Xcl is the mass of Xci and e is the xa polarization vector. The decay constant f Xcl is 
a nonperturbative quantity and may be estimated from the potential models, the QCD sum 
rules, or lattice QCD calculations. 

The leading-twist light-cone distribution amplitude of Xci can be expressed as 

(Xci( P ,e)\c a (y)cp(z)\0) = £ jfW e 4 '^ 1 ^^ (5) 

where u and 1 — u are respectively the momentum fractions of the c and c quarks inside the 
Xci meson, and the wave function Xcl (u) for x& meson is symmetric under u <-> 1 — u. 

In naive factorization, we neglect the strong interaction corrections and the power correc- 
tions in Aq CD /7?v Then the decay amplitude at leading order is written as 

iM = 2if Xcl m Xcl e* ■ PB F x {mlJ^ [WE(tf a + ^) - V tb V t l(C 3 + ^ - C 5 - ^)], (6) 
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Figure 1: Feynman diagrams for vertex and spectator corrections to B — > x&K- 

where ps is the momentum of B meson, F\ is the B — > K transition form factor and N c is 
the number of colors. We do not include the effects of the electroweak penguin operators since 
they are numerically small [TT] . 

The form factors for B —* K are given as 



(K(Pk) \s^b\B(p B )) = (p B + p K )^ - 



i i 
m B — m K 

p2 



Flip 



2 „ , m\-m\ 
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P»F (p 2 



(7) 



where p = p B — Pk is the momentum of Xci with p 2 = m Xcl , m B and uik are respectively the 
masses of B and K mesons. We will neglect the kaon mass for simplicity. We can use the ratio 
between these two form factors as 



F (p 2 
F^p 2 



1 



P 



m 



2 ' 

B 



So we need only one of the two form factors to describe the decay amplitude. 

Taking nonfactorizable corrections in FigHinto account, the full amplitude for B — > x c \K 
in QCD factorization is written compactly as 



Gp 



%M = 2if Xcl m Xcl e* ■ p B F l (m Xcl )-y=[V cb V*a 2 - V tb V t * s (a 3 - a 5 



(9) 



where the coefficients (i = 2, 3, 5) in the naive dimension regularization(NDR) scheme are 
given by 



• c ^ " 1 ' c 



as = C7 8 + ^ + ^C7 4 (-18 + 121n^ + / J + / 

m b 



u 



(10) 



C 6 a s C F / ni b 
a 5 = C 5 + ----C 6 (-6 + 121n- + // + /j 



n 



where C F = (N 2 — 1)/(2N C ) and m b is the b quark mass. 
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The function // is calculated from the four vertex correction diagrams (a, b, c, d) in Fig^ 
and reads 



r± roll Zu 

fi= du(f) x (u) lnw — 3«7r + 3 ln(l — z) + 

jo c L 1 — u 



m 



(1-zu) 2 



+— -+b rr x)z 2 uln(zu , 11 

where z = m Xci /m B , and we have already symmetrized the result with respect to u «-> 1 — w. 

The function /// is calculated from the two spectator correction diagrams(e, f) in FigHJand 
it is given by 

/// = 2 / d£— t— / <fy / du(f) Xcl {u R- + -7Z r, 12 

N c Fi(m Xol )m B 1 — z Jo £ Jo y Jo u y{l — z) 

where £ is the momentum fraction of the spectator quark in the B meson and y is the momentum 
fraction of the s quark inside the K meson, <f) B , <pK are the light-cone wave functions for the B 
and K meson respectively, fn corresponds to the T 11 term in Eq.flTJ. 

The spectator contribution depends on the wave function 0# through the integral 

'^MU^. (13) 

Since 0_b(£) is appreciable only for £ of order Aqcd/wib, is of order Aqcd- We will choose 
\b ~ 300 MeV in the numerical calculations [3]. 

There is an integral related to <pK{y) in Eq. (jl2j) which will give logarithmic divergence. 
Therefore QCD factorization breaks down even at leading order. This is different from B — > 
J/ipK decay which does not have logarithmic divergence at leading twist [3 |Bj. The reason is 
that the logarithmic divergences arising from contributions of vector and tensor currents are 
cancelled out in the B — > J/ipK decay, whereas there is no such cancellation for the B — > Xc\K 
decay. Following Ref . , we treat the divergent integral as an unknown parameter and write 

\t^ = t\«^zA = mH . 1)t (14) 

y z Jo y 

where the asymptotic form (pxiy) = %(1 — y) is used for the kaon twist-2 light-cone distribu- 
tion amplitude. And we will also choose X H = \n.{rriBl h-QCD)® as a rough estimate in our 
calculation. 

For numerical analysis, we choose i<\(m 2 cl ) = 0.743^3] and use the following input param- 
eters: 

m b = 4.8 GeV, m B = 5.28 GeV, m Xcl = 3.5 GeV, 
f Xcl = 230 MeV, f B = 180 MeV, f K = 160 MeV. (15) 

Here for the Xci decay constant f Xcl we give its expression in the heavy quarkonium potential 
model as 



V ™>Xcl 

where R'(0) is the derivative of the Xci radial wave function at the origin, and we may roughly 
estimate its value by using some potential models, e.g. given in Ref.|14j. 
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c a 


c 2 


c 3 


c 4 


c 5 


c 6 


NDR 


1.082 


-0.185 


0.014 


-0.035 


0.009 


-0.041 



Table 1: Next-to-leading-order Wilson coefficients in NDR scheme (See Ref.jTTj) with \x = 4.4 
GeV and A^L = 225 MeV. 





(«) 


a 2 




a 3 


a 5 


6u(l 


-«) 


0.1255-0.0^ 


}15i 


0.0040+0.0026i 


-0.0027-0.0031i 


S(u- 


1/2) 


0.1154-0.0^ 


}14i 


0.0043+0.0026i 


-0.0031-0.0031i 



Table 2: The coefficients tij at [A — 4.4 GeV with different choices of 4> Xcl {u). 



The specific form of the Xci light-cone distribution amplitude is not known, but as in the 
J/ip case (see |7J |H|) we may use the light meson asymptotic form 4> Xcl {u) = 6u(l — u) as 
an approximate description. In the numerical analysis, we also consider the form <f) Xcl {u) = 
S(u — 1/2), which comes from the naive expectation corresponding to the nonrelativistic limit 
of the distribution amplitude. Although there are uncertainties associated with the form of the 
distribution function, we will see shortly that the calculated decay rates are not sensitive to 
the choice of the distribution amplitude. 

Using the next-to-leading-order (NLO) Wilson coefficients in Table. ^ we get the results of 
coefficients cij which are listed in Table. El With the help of these coefficients a«, we calculate 
the decay branching ratios with two different choices of the Xci distribution amplitude and get 

Br(5^ Xcl K) = 0.97xl0" 4 , (0 Xcl («) = 6u(l - «)); (17) 
Bi(B -> XdK) = 0.87 x 10- 4 , (<f> Xcl (u) = 5(u - 1/2)). (18) 

The measured branching ratio is[9j 

Br(B° -> XdK ) = (5.4 ± 1.4) x 10" 4 , (19) 

which is about six times larger than the theoretical results. 

We now consider the B — > XcoK decay. Unlike the x& mode that we just analyzed, Xco does 
not have couplings to V or A currents, so this decay mode is prohibited in naive factorization. 
But it may occur if there is an exchange of an additional gluon^HJ El- The branching ratio 
from experiment isfTU] 

Br(B + -> XcoK + ) = (6.0 ± 2.1) x 10" 4 . (20) 

Similar to the calculation performed above, we write the Xco light-cone distribution ampli- 
tude in the following general form 

(Xco( P )\c a (y)cp(z)\0) = \ [ due^y + ^- u ^ [f Xc0 m Xc0 l^ x Ju) + flM P Jl»] > (21) 



o 



where f Xc0 ,f Xc0 , ( fi Xc0 ( u ) an d 4> Xc0 ( u ) denote the Xco decay constants and light-cone distribution 
amplitudes, respectively for the scalar and vector currents. It is easy to see that for the 
scalar meson Xco the decay constants for the vector current and the antisymmetric tensor 
current all vanish (i.e., (Xco|c7^c|0) = (Xcolccr^clO) = 0), and only the scalar current decay 
constant is nonzero. Therefore on the right hand side of Eq. (J21j) only the first term makes a 
real contribution and it is essentially the leading twist (twist-3) contribution. However, here 
we also list the second term (twist-2) in Eq. (j21J) in order to show that it would also give an 
infrared divergence contribution if the vector current decay constant were nonzero. 
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Because of charge conjugation invariance, 4> s Xc0 {u) is symmetric while <p XcQ (u) is anti- 
symmetric under u <-> 1 — u. The Feynman diagrams of order a s correction for B — > XcoK 
are the same as those in FigO In the calculation of Xci, the contribution of the four vertex 
diagrams in Fig. ^ is infrared safe. However, we find there are infrared divergences arising from 
the vertex diagrams in the B — > XcoK decay. Keeping only the divergent terms, we get 



x [- hi A In z + (1 - z)[\n A + i ln 2 (^) - m H^)]} }, (22) 



where z = m Xc0 /m% and A corresponds to the soft gluon momentum cutoff divided by the 
charm quark mass. The infrared divergence is explicitly seen as A — > 0. 
The contribution arising from spectator interactions is given by 

<x,C F ir. t f i m c m Xc0 fS f 1 ^ c <Pb{0 f 1 ^.Mv) f 1 ^, <^co( M ) 



l m B — m x J o £ Jo y z Jo u(l — u) 



+r M f\*sm l l dy tM. f du fkM]. (23) 

Jo 4 Jo y Jo 1 - u J 

There are logarithmic divergences even in the leading-twist order. So, soft gluon exchange 
dominates B — > XcoK decay. 

In the above calculations, we have chosen the heavy quark limit as mj, — > oo with m c /mb 
fixed. In this limit, the logarithmic divergences in Eq. ()12j1 for B — > Xc\K decay, and the 
infrared divergences in Eq.()22|) as well as the logarithmic divergences in Eq.()23p for B — > XcoK 
decay, will break down QCD factorization. Another choice of the heavy quark limit is that 
rrib — > oo with m c fixed. Then all the divergences mentioned above are power corrections and 
should be dropped out, so QCD factorization still holds in this limit. Physically, the latter 
case is equivalent to the limit of zero charm quark mass in which charmonium is regarded as 
a light meson. Obviously, the first choice of the heavy quark limit with m c /m^ fixed is more 
relevant in phenomenological analyses, and is also usually used in theoretical studies [3 HI] , 
where it is expected that QCD factorization should apply to B meson exclusive decays into 
charmonium in the limit m c — > oo with corrections of order Aqcd/(^c«s) ~ 1- Our result 
shows that this expectation holds for decays like B — > J/ipK but not for e.g. B — > XcoK where 
the vertex infrared divergence will break down factorization at order of A.qcd/ (m c a s ). It is 
also worthwhile to emphasize that the infrared divergences in Eq.()22j) are more serious than 
the logarithmic divergences in Eqs. (fT2^) . (f2*3^ . This is because, when one includes the parton 
transverse degrees of freedom, end-point singularities arising from spectator interactions can be 
regularized within the framework of kx factorization ^H] and logarithmic divergences can then 
be removed. However, the infrared divergences due to nonfactorizable vertex corrections still 
exist, and they are independent of the specific form of the light-cone distribution functions. 

In summary, we have studied the exclusive decays of B meson into P-wave charmonium 
states Xcj{J — 0, 1) in the QCD factorization approach with light-cone distribution functions 
describing mesons involved in the processes. For B — > x&K decay, the factorization breaks 
down due to logarithmic divergences arising from nonfactorizable spectator interactions at twist- 
2 order and the decay rate is too small to accommodate the data. The situation for B — > XcoK 
decay is even more serious because of the infrared divergence arising from nonfactorizable vertex 
corrections as well as logarithmic divergences due to spectator interactions at the leading-twist 
order. Considering the above problems encountered in B — >• XcjK decays and especially the 
infrared divergence in B — > XcoK decay, we intend to emphsize that the QCD factorization 
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method with the present form may not be safely applied to exclusive decays of B meson into 
charmonium, e.g. P-wave states. However, it should be noted that our results are obtained 
in QCD factorization by using the light-cone distribution functions for charmonium (as in 
e.g. refs. j7[ El), which may not be a very appropriate description for heavy quarkonium. To 
further clarify the problem, many theoretical studies such as use of the nonrelativistic QCD for 
charmonium, the perturbative QCD factorization method, and the nonperturbative QCD effects 
should be attempted. In any case, a better understanding is definitely needed for describing 
.B-meson exclusive decays to charmonium, and for reducing the big gap between the present 
theoretical results and the experimental observations. 
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